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ABSTRACT 

Context. A large fraction of stars are found in binary systems. It is therefore important for our understanding of the star 
formation process, to investigate the fragmentation of dense molecular cores. 

Aims. We study the influence of the magnetic field, ideally coupled to the gas, on the fragmentation in multiple systems of 
collapsing cores. 

Methods. We present high resolution numerical simulations performed with the RAMSES MHD code starting with a uniform 
sphere in solid body rotation and a uniform magnetic field parallel to the rotation axis. We pay particular attention to the 
strength of the magnetic field and interpret the results using the analysis presented in a companion paper. 
Results. The results depend much on the amplitude, A, of the perturbations seeded initially. For a low amplitude, A = 0.1, we 
find that for values of the mass-to-flux over critical mass-to-flux ratio, p, as high as n — 20, the centrifugally supported disk 
which fragments in the hydrodynamical case, is stabilized and remains axisymmetric. Detailed investigations reveals that this 
is due to the rapid growth of the toroidal magnetic field induced by the differential motions within the disk. For values of p 
smaller ~ 5, corresponding to larger magnetic intensities, there is no centrifugally supported disk because of magnetic braking. 
When the amplitude of the perturbation is equal to A = 0.5, each initial peak develops independently and the core fragments 
for a large range of [i. Only for values of jj, close to 1 is the magnetic field able to prevent the fragmentation. 
Conclusions. Since a large fraction of stars are binaries, the results of low magnetic intensities preventing the fragmentation in 
. £h ' case °^ wea, k perturbations, is problematic. We discuss three possible mechanisms which could lead to the formation of binary 
systems namely the presence of large amplitude fluctuations in the core initially, the ambipolar diffusion and the fragmentation 
» , ■ during the second collapse. 

d; 

Key words. Magnetohydrodynamics - Instabilities - Interstellar medium: kinematics and dynamics - structure - clouds 



1. Introduction 

Understanding the fragmentation of collapsing prestellar 
dense cores, is of great importance to our understanding of 
the star formation process. In particular, determining the 
number of fragments, their masses and their orbital char- 
acteristics are fundamental and rather challenging prob- 
lems. Most of stars are found in binary or multiple sys- 
tems (Duquenoy & Mayor 1991). Many studies have in- 
vestigated this issue in the context of hydrodynamical 
calculations (e.g. Miyama 1992, Boss 1993, Bonnell 1994, 
Truelove et al. 1998, Bate & Burkert 1997, Bodcnheimer 
et al. 2000, Matsumoto & Hanawa 2003, Hennebelle et al. 
2004, Goodwin et al. 2004, Banerjee et al. 2004). Although 
detailed conclusions appear to be very sensitive to the ini- 
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tial conditions (i.e. core shape, thermal, rotational and 
turbulent energy as well as to the equations of state), 
some trends can nevertheless be inferred. It is widely ac- 
cepted that, under realistic initial conditions, a collapsing 
dense core fragments into a few objects, the exact number 
depending on the specific conditions. Therefore, fragmen- 
tation of a rotating collapsing dense core appears to be 
the most widely accepted mechanisms to explain the for- 
mation of multiple systems. Indeed, it is today the only 
viable mechanism since other possibilities such as fission 
of a rotating protostar or capture of a companion, fail in 
realistic conditions to produce a large fraction of binaries 
(Bodenheimer et al. 2000). 

Due to the difficulty of the problem and despite its im- 
portance, the question of the role of the magnetic field in 
this process has remained little addressed. With the recent 
progress achieved in numerical techniques as well as the 
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increasing computing power, several studies have recently 
investigated this issue. Hosking & Whitworth (2003) us- 
ing an SPH two fluids code, conclude that magnetically 
subcritical cores do not fragment. Machida et al. (2005) 
using a nested grid code performed an extensive num- 
ber of calculations varying the initial core rotation and 
magnetic field. They find that fragmentation is possible 
if the rotation is sufficiently large and the magnetic field 
strength sufficiently small. Fragmentation of a collapsing 
magnetized cloud using adaptive mesh refinement tech- 
niques have also been studied by Ziegler (2005), Bancrjcc 
& Pudritz (2006) and Fromang et al. (2006) which all find 
that the magnetic field has a strong influence. Price & 
Bate (2007), using magnetized SPH techniques which in- 
sure the nullity of the divergence of the magnetic field, find 
that with large initial perturbations, even for large values 
of the magnetic strength, the fragmentation is possible. In 
all these studies, it has been found that the magnetic field 
has a strong impact on the fragmentation of the collapsing 
dense core. 

In this paper, we study further the fragmentation of 
collapsing cores, focusing mainly on the influence of the 
magnetic field strength and the amplitude of the initial 
perturbation. Using the analysis developed in the compan- 
ion paper by Hennebelle & Fromang (2007), thereafter pa- 
per I, we try to understand, when possible quantitatively, 
the physical reasons for the various numerical results. As 
in paper I, we consider initially a uniform density cloud 
in solid body rotation, threaded by a uniform magnetic 
field. At this stage, we restrict the problem to the case 
where magnetic field and rotation axis are aligned. We 
note that Price & Bate (2007) have investigated the in- 
fluence of initially perpendicular rotation axis and mag- 
netic field. More realistic initial conditions, including non 
uniform density field, turbulent velocity field and rotation 
axis non aligned with the magnetic field, will be considered 
in future studies. In particular, we determine the lowest 
value of /i, the mass-to-flux over critical mass-to-flux ratio 
for which fragmentation is suppressed. It turns out that 
the results depend strongly on the amplitude of the initial 
perturbation. 

The paper is organized as follow. In the first section, 
we consider the same initial conditions as used in paper I, 
namely a spherical uniform one solar mass dense core. The 
thermal over gravitational energy of the core is equal to 
~0.37 and the rotation over gravitational energy is 0.045. 
Note that such rotation are typical of (may be slightly 
larger than) values observed in dense cores (Goodman et 
al. 1993). The initial density is about ~ 5 x 10~ 18 g cm~ 3 
and the cloud radius, i?o —0.016 pc. The freefall time is 
thus of the order of 3 x 10 4 years. We use a barotropic 
equation of state: C 2 S = (C°) 2 x (1 + (p/ p^ 4 / 3 ) 1 / 2 , where 
C s ~ 0.2 km/s is the sound speed and p c = 10~ 13 g cm~ 3 . 
This set of cloud parameters is known to give rise to a 
disk that unambiguously fragments in the pure hydrody- 
namics case. The only difference with paper I is that an 
r7i = 2 perturbation of amplitude A = 0.1 in the density 
field, p(r,6,z) = po x (1 + Acos(m6)), as well as in the 
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Fig. 1. Equatorial density and velocity field for mass-to- 
flux over critical mass-to-flux ratio p — 1000 and pertur- 
bation amplitude equal to A = 0.1. The size of the figure 
corresponds to about one tens of the cloud initial radius. 
In physical units, this is about 300 AU. The freefall time 
is about 3 x 10 4 years. 

B z field, is added. In the second section, we further dis- 
cuss the physical interpretation of the results obtained in 
our simulations and we estimate analytically the critical 
value of p for which it is expected that the magnetic field 
stabilizes the disk. In the third section, we consider ini- 
tial density perturbations of amplitude 0.5 and show that 
in this case, fragmentation can be obtained for a much 
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Fig. 2. Same as Fig. Q]for /x = 50 and A = 0.1. 



Fig. 3. Same as Fig. Q] for // = 20 and A = 0.1. 



wider range of /i. The fourth section provides a discussion 
on possible mechanisms leading to disk fragmentation and 
the formation of binaries, even in the presence of magnetic 
fields. The sixth section concludes the paper. 

2. Weak initial perturbations 

We investigate the fragmentation in the case where the 
amplitude of the perturbations is initially A — 0.1. In this 
case it is found that without rotation no fragmentation 
occurs. Therefore, we call this type of initial condition, ro- 
tationally driven fragmentation. We next consider various 



values of the magnetic mass-to-flux over critical mass-to- 
flux ratio, namely \i = 1000 (quasi hydrodynamical case), 
H = 50, 20, 5, 2 and 1.25. Recalling that fi = 1 corre- 
sponds to the case of a cloud supported by the magnetic 
field, the last value corresponds to a strongly magnetized 
supercritical cloud. 

In the following, we display the inner part of the col- 
lapsing cloud. The size of each plot is about ~ i?o/10. This 
corresponds to a physical size of 1.6 x 10 -3 pc or about 
-300 AU. 
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Fig. 4. Same as Fig. [T] for p, = 5 and A = 0.1. 



2.2. Hydrodynamic case: fragmentation 

We present the hydrodynamical case (/U = 1000). Figure Q] 
shows three snapshots. Equatorial density and velocity 
field are displayed. The first snapshot is taken shortly after 
the formation of the first Larson core. A rotationally sup- 
ported structure has formed and since it is very unstable, 
a strong spiral pattern is developing. This is very similar 
to the results obtained by other authors (e.g. Matsumoto 
& Hanawa 2003, Hennebelle et al. 2004, Goodwin et al. 
2004). The next snapshot shows that the central bar-like 
structure, fragments into two objects whereas the spiral 
pattern keeps growing and expanding due to further ac- 
cretion of mass and momentum. The last snapshot shows 
that the centrifugally supported structure continues frag- 
menting and that 5 fragments have developed. Altogether, 
this is very similar to the results reported by many authors 
investigating the fragmentation of hydrodynamical rotat- 
ing cores. We note that in spite of the strongly non-linear 
evolution, the result remains remarkably symmetric due to 
proper resolution of the Jeans length (10 numerical cells 
per Jeans length). 



2.2. Weak magnetic field cases: suppression of 
fragmentation 

We present results for the weak field cases, p = 50 and 
p = 20. Figure [2] shows results for p = 50. Due to the dif- 
ference in the initial magnetic strengths, it is not possible 
to compare the simulations at exactly the same times. The 
first snapshot shows a spiral pattern which is similar to the 
one seen in the previous case. A difference however, is that 
the central bar-like structure appears to be much shorter. 
This is likely due to the central magnetic field which is 
strongly amplified by the rapid twisting of field lines. As 
a result, the cloud first fragments into three fragments in- 
stead of two. This is due to the fact that since the central 
bar-like structure does not fragment in two objects as in 
the previous case, more material and angular momentum 
is available in its vicinity to produce one object on each 
side. The next snapshots show that a symmetry breaking 
occurs due to one of the two satellites having merged with 
the central more massive fragment. The structure in the 
outer part (x > 0.03i?o) is still quite symmetric. Further 
fragments are forming in the outer part of the spiral pat- 
tern at x = ±0.05i?o- 

From these results, one concludes that even for values 
of p as large as 50, the magnetic field has a significant im- 
pact on the evolution of the centrifugally supported inner 
structure particularly its fragmentation. This is due to the 
strong amplification of the toroidal and radial magnetic 
field generated by the differential motions in the collaps- 
ing core and in the centrifugally supported structure (see 
paper I and following sections). 

Figure [3] shows three snapshots for p = 20. The first 
snapshot shows a weak spiral pattern in the inner part. 
The second and third panels shows that the centrifugally 
supported structure grows as in the previous cases but 
it remains much more uniform and the spiral pattern is 
much less pronounced. As shown in paper I, the angu- 
lar momentum profile is similar to the hydrodynamical 
case since the magnetic braking is very weak during the 
collapse phase. Therefore, the differences with the case 
p = 50, is not smaller angular momentum but rather a 
stronger magnetic field. Indeed, with a magnetic toroidal 
component, the velocity of the fast MHD wave, which can 
be loosely seen as a effective sound speed, is: \/C% + 
where V Q 2 = Bg/Airp (see section [3]). Therefore, the disk 
is stabilized against gravitational fragmentation. As a re- 
sult, no fragmentation is obtained, only one central star 
forms and grows by accretion. We stress that the value of 
the magnetic field corresponding to p — 20, is very mod- 
est, far below the values inferred from observations which 
indicates that \i = 1 — 2 (Crutcher 1999). 

Finally, to explore whether the suppression of the frag- 
mentation depends on initial conditions, such as the values 
of a and (3 (the ratios of thermal and rotational over grav- 
itational energies), we have performed runs with different 
a and (3 values. First, we explore the influence of stronger 
and weaker rotation choosing (3 = 0.2 and (3 = 0.02 re- 
spectively keeping the value of a constant. The behaviour 
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of these runs is very similar to the one with (3 = 0.045 
presented here. A large disk forms but it remains well ax- 
isymmetric and does not fragment. We have also explored 
the effect of a smaller thermal energy, taking a = 0.2 and 
keeping (3 = 0.045. The disk is a little more axisymmet- 
ric than in the case a = 0.37 but still does not fragment 
either. 

We conclude that, for small initial density perturba- 
tions, a magnetized dense core ideally coupled to the mag- 
netic field and having /i = 20 does not fragment for a large 
range of initial conditions. 

It is worth to compare with the study of Machida et 
al. (2005) although it is not straightforward because their 
initial conditions consist in a filament and therefore are 
different from ours. Also they used different definitions 
to quantify the amount of rotation and magnetic ener- 
gies in their simulations. Comparison with their Fig. 10 is 
worthwhile nevertheless. Their parameter lo turns out to 
be equal to \fj3. However, their initial state has a peak 
density of about 10 4 g cm -3 which is less dense than 
ours by a factor of roughly 100. Since the ratio of ro- 
tational over gravitational energies increases during the 
collapse, we can say that our (3 = 0.045 correspond to 
uj < 0.2 in Machida et al. study (for a homologous con- 
traction (3 oc 1/r cx /o -1 / 3 , this would indicate that (3 
should be divided by roughly ~ 4 — 5 to be compared 
with Machida et al.'s results). For values of co < 0.2, they 
found that for B z / \J%i:C1p c > 0.1, they have no fragmen- 
tation. Comparing this value with our parameters is again 
not straightforward but B z / \f%irUfp~ c is about ~ 1/ 'y/afi. 
For a = 0.37 and \i = 20, l/y/a/j, ~ 0.1. Therefore, our 
results broadly agrees with the results of Machida et al. 
(2005). 

Comparison with the study of Price & Bate (2007) is 
not possible at this stage, since they use a barotropic equa- 
tion of state which becomes adiabatic at 10 -14 g cm' 3 . 
Thus, even the hydrodynamical calculations they present 
does not fragment when the perturbation is weak (their 
figure 3). 

2.3. Intermediate magnetic intensity: formation of 
pseudo-disks 

We present results for the intermediate magnetized cases, 
namely fi = 5 and fi — 2. Figure [5] shows two snapshots 
for /i = 5. As shown in paper I, no centrifugally supported 
structure forms instead a magnetized pseudo-disk devel- 
ops. Pseudo-disks arise when disk-like structures form, 
say oblate ellipsoids, which are not supported by rotation 
but rather by magnetic support. According to the analy- 
sis presented in paper I, the angular momentum is lower 
for these values of fi primarily because the collapse occurs 
first along the field lines. Thus, the material within the 
pseudo-disk and central object was initially located along 
the pole and has less angular momentum. Some magnetic 
braking also occurs, reducing the angular momentum fur- 



ther. As a consequence, no centrifugally supported disk is 
observed and no fragmentation is occurring. 

Figure [5] shows two snapshots for [i = 2. The first 
panel shows that a filamentary structure has developed 
after t = 1.51 freefall times. We believe that this filament 
is due to the non-linear evolution of the initial to = 2 per- 
turbation. The reason why it has a shape different than 
in the previous case is that since the cloud is more sup- 
ported by the magnetic field, the collapse lasts 30% longer. 
Therefore the perturbation has more time to develop and 
to become non-linear. Indeed, if no perturbation is initially 
included, an axisymmetric pseudo-disk develops. 

2.4. Nearly critical core 

We now consider the case of a very magnetized supercriti- 
cal core having fj, = 1.25. As expected the collapsing time 
is now longer and roughly equal to two freefall times. As 
revealed by the two panels of Fig. [6l the collapsing dense 
core remains almost axisymmetric in spite of the initial 
to = 2 density perturbations. This is due to the strong 
magnetic support which prevents the development of the 
perturbation. Indeed, since the cloud is nearly supported 
by the magnetic field, the magnetic Jeans Mass is just 
slightly smaller than the cloud mass. Therefore, each peak 
of the to = 2 density perturbation is gravitationally stable. 
Moreover, unlike the thermal support, the magnetic sup- 
port does not decrease during the collapse since the ratio 
of magnetic over gravitational energy stays roughly con- 
stant. Therefore the magnetic Jeans mass stays roughly 
constant, unlike the thermal Jeans mass. 

3. Physical interpretation 

One of the important conclusions reached in the previous 
section is that fragmentation can be suppressed even for 
large values of fj,, i.e. weak magnetic fields. Here we dis- 
cuss the physical reason for the disk stabilization arising 
at low magnetic strength. We stress that since the angular 
momentum is nearly identical for /i = 20 to the hydrody- 
namical case (see paper I), the reason of this stabilization 
is not magnetic braking. 

3.1. Influence of the magnetic field on the disk stability 

The influence of the magnetic field on the stability of 
a self-gravitating disk, has first been investigated by 
Lynden-Bell (1966) for a uniform rotation. The dispersion 
relation he obtained (see his Eq. 1) entails an effective 
sound speed \/Cf+V^ showing that the magnetic pres- 
sure term has a stabilizing influence. However, there is also 
a destabilizing contribution due to the magnetic tension 
term. As a result, such configurations are unstable. 

Elmegreen (1987) and Gammie (1996) both consider 
the influence of the magnetic field on the stability of a dif- 
ferentially rotating system. In that case, the shear drives 
the growth of a toroidal magnetic component which sta- 
bilizes the disk. They both conclude that whereas in the 
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Fig. 5. Same as Fig. [T] for /z = 2 and A = 0.1. 



Fig. 6. Same as Fig. Q] for p = 1.25 and A = 0.1. 



absence of shear the magnetic field is strongly destabiliz- 
ing, it has a stabilizing influence when significant shear 
is present. More precisely, Gammie (1996) computes the 
response of the disk to nonaxisymmetric perturbations for 
various value of Q, the Toomre parameter. He shows that 
the response is much weaker when a substantial magnetic 
field is present and concludes that the stabilizing effect of 
the field becomes significant once the magnetic pressure 
is comparable to the gas pressure. This is qualitatively in 
good agreement with our results. In particular, we see that 
the disk in case /z = 20 remains much more axisymmetric 
than in the case p = 1000 and p = 50. 

In order to quantitatively verify that this mechanism is 
responsible of the disk stabilization, we have computed the 
Alfven speed associated with the azimuthal component of 
the magnetic field, Bg/ ' ^JAirp, within the disk at various 
times corresponding to the period during which fragmen- 
tation is taking place. Since both Bq and p vary, we have 
integrated at every radius along the z-axis through the 
disk. As explained in paper I, as soon as the azimuthal 
component of the magnetic field exceeds the ram pressure 
of infalling material, a magnetic tower builds up and the 
disk starts expanding. In this case, the disk boundaries are 
not clearly defined. We adopt here a density threshold to 
define the disk and stop the integration when p < 10~ 14 
g cm -3 . We calculate < >— J Bg/Anp x pdzj J pdz, 



which allows us to estimate + V^/C^ within the disk. 
Figure [7] shows y/l + Vf/C* for p = 1000, 50 and 20. 
As can be seen, for the case p = 1000, the quantity 
•v/l+ < Va > /Cg remains smaller than 1.2 at all the 
times displayed. Figure. Q] shows that by this time, frag- 
mentation has already occurred. The second and third 
panel of Fig. [7] show that < V£ > /Cf grows much 

more rapidly for smaller p reaching values greater than 2. 
The growth is more rapid for p = 20 than it is for p = 50. 
We conclude that the growth of the toroidal component 
induced by the differential rotation, is mainly responsible, 
for disk stabilization. 

For simplicity, the works mentioned above, have re- 
stricted the analysis to a thin disk. However, the expan- 
sion of the magnetic tower triggered by the growth of the 
toroidal magnetic field, removes some material from the 
disk and therefore reduces the disk surface density. This 
makes the disk even less prone to fragment. This effect 
which is not taken into account in the thin disk analy- 
sis, certainly contributes to further stabilize the centrifu- 
gally supported structure against fragmentation. Indeed, 
the mass of the disk and the mass within the tower, turn 
out to be roughly comparable. 
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Fig. 7. v/l + KVCI within the disk. 



3.2. Analytical estimate of the critical value of \i 

We estimate analytically the value of [i at which one ex- 
pects to find a significant influence of the magnetic field on 
the fragmentation. Following Gammie (1996), we will con- 
sider that a strong stabilizing influence is achieved when 
the Alfven speed associated to the toroidal component is 
comparable to the sound speed. 

The principle of the analysis is as follow: we compute 
the growth rate of the toroidal component inside the cen- 
trifugally supported structure, so that we can estimate the 
time, T mag , needed for the Alfven speed to become compa- 
rable to the sound speed. We then compute the dynamical 
time of the disk, Td yn , over which fragmentation will oc- 
cur. The critical value of fx, below which fragmentation is 
quenched, is obtained when Td yn is equal to r mag . 

The growth time of the toroidal magnetic component 
within the disk has two contributions. First the twisting of 
the radial component due to differential rotation propor- 
tional to VeB r /r and second the wrapping of the vertical 
magnetic component due to the vertical gradient of B z 
proportional to VgB z /z, where z is the disk height. The 
divergence constraint shows that B r /r ~ B z /z, indicating 
that the 2 contributions are comparable. 

Thus, the growth of Bg can be estimated by writing 
d t Bg ~ B z Vg/h where h is the disk height. Since we are 



interested in computing the growth of the magnetic field 
from initially small values, we neglect the influence of the 
magnetic tower on the disk vertical expansion and write: 
h ~ ^C 2 J{AirGp). With Eq. (1) of paper I, we get h ~ 
r/V2d, where d is the ratio of the density over singular 
isothermal density. 

If we assume that the disk is roughly Keplerian, we 
get Vg ~ y/M s G/r, M s being the mass of the star and 
the disk. Taking the expression of B z stated by Eqs. (1) 
of paper I, we obtain: 



Bg~tx 
The 



r 5/2 

value of 



(1) 

is obtained by requiring 



(Bg I \f^Kp) jC a ~ 1 . Using the expression of density stated 
by Eq. (1) of paper I, we obtain: 

r 3/2 

Tmag " H Z VGW S - (2) 

The dynamical time of the disk, i.e. the time relevant 
for fragmentation, is the rotation time: 



27T7- 3 / 2 



Tdyn 



(3) 



(4) 



Vgw s ' 

Therefore, we obtain: 

Tmag ^ ! 
Tdyn 2nH z 

To compute the value of H Zl one can use Table 1 of pa- 
per I, but it is more convenient to express the result in 
terms of the initial cloud parameters namely /i and a. 
The latter is the ratio of thermal over gravitational en- 
ergy and is given by a — 5/2Cf Rq/(MG). From Eq. 4 
of paper I, the value of H z is given by VGB^R /(2Cg), 
whereas from the definition of /i, B® can be written 
as M /(irRl)^^^^/{Q.b?>ii). Gathering the different 
expressions, we get: 

Tmae "\/2 . . 

— * - irjT x m ' (5) 

Tdyn 3V57T 

We can deduce the critical value of fi 
15 

Mc - — 



(6) 



Note that this value is approximate and should be valid 
within a factor of a few. 

For a ~ 0.37, we obtain [i c ~ 40. This analytical esti- 
mate is therefore in reasonable agreement with our nu- 
merical simulation, since for /i — 20, fragmentation is 
suppressed. In this case, the toroidal magnetic pressure 
becomes comparable to the sound speed in about half a 
rotation period. Since typically fragmentation is occurring 
over a few rotation period, such a fast growth time appears 
to be sufficient to prevent the disk from fragmenting. On 
the other hand for /i = 50 the disk is fragmenting, con- 
firming that this case is slightly above the critical value. 

For n = 1000, one has T mag /Td yn — 20. Therefore in 
such a case, the growth of the toroidal component is too 
slow to significantly influence the disk evolution. 
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Fig. 8. Same as Fig. Q] for \i = 20 and A = 0.5. Note nev- 
ertheless that the figure corresponds to a length of ~1200 
AU, four times bigger than the length displayed in Fig. Q] 



4. Strong initial perturbations 

In this section, we investigate the effect of a density fluc- 
tuation having a larger initial amplitude, A = 0.5. In that 
case the cloud is more prone to fragment because the per- 
turbations develop quickly. Such strong perturbations can 
be due either to initial conditions or to strong external 
triggering. Although we are keeping in this section the 
same amount of rotation as in the previous section, we 
also performed simulations with the same initial condi- 
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Fig. 9. Same as Fig. [T] for fj, = 2 and A = 0.5. 



tions but with no rotation and found similar results during 
the early phase of fragmentation as in the case with rota- 
tion. This indicates that with such strong perturbations, 
the initial cloud fragmentation is purely thermal and in- 
dependent of the rotation. We refer to this situation as 
thermal fragmentation. 

4.1. Thermal fragmentation 

We present results for weak and intermediate values of 
magnetic strengths, namely /i — 20 and \i = 2. Figure [8] 
shows results for [i = 20. The first panel shows that a fil- 
amentary structure develops in which two protostar have 
formed. The two protostars are well separated by roughly 
~ 0.15 x Rq as opposed to ~ 0.01 x R in the case of weaker 
perturbations. As explained previously, the fragmentation 
in this case has a purely thermal origin and is simply due 
to the development of the initial perturbations (in the case 
of weak perturbations, rotationally driven fragmentation 
docs not strongly depend on the perturbation amplitude). 
The second panel shows that the two fragments are ap- 
proaching each other falling along the filament. The final 
panel shows that the two fragments have merged and that 
a protostellar disk forms. At this point, however, it is nec- 
essary to remember that the second collapse is not treated 
in the simulation and that only the first Larson core is con- 
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Fig. 10. Same as Fig. Q] for fx = 1.25 and A = 0.5. 

sidercd. Therefore the merging should be considered with 
great care. Careful thermodynamical treatment of the first 
Larson core is required. Indeed the first core is known to 
have a short life time of roughly 10 3 years. Therefore since 
the timeshift between the first and second panel is larger 
than 10 4 years, it is expected that second collapse would 
have already occurred by the time, the two fragments ap- 
proach each other. Since stars have very small radii, it 
is unlikely that the binary is going to merge. Therefore 
the final product could be a binary rather than a single 
protostar. 

Figure [5] shows results for n — 2. Both panels show 
that two fragments develop in a similar way as for fi = 20. 
Unfortunately, the relative strong value of the magnetic 
intensity, makes the time steps much shorter than in the 
previous case and it becomes computationally expensive 
to follow the long range evolution. It seems however likely 
that a similar evolution is expected leading to the merging 
and the two first cores. 

4.2. Strong field: suppression of thermal fragmentation 

Figure [TU] shows results for /i = 1.25 and A = 0.5. As 
revealed by the first panel, the m = 2 density perturbation 
has only weakly developed and only two shallow density 
maximum are visible. The second panel shows a single 



object in the centre surrounded by a pseudo-disk. Since 
the timeshift between the two panels is about 0.015 freefall 
time (~ 5 x 10 2 years) there is not sufficient time for the 
second collapse to arise. It seems plausible that in this 
particular case, a single star is likely to form. 

Therefore, it appears, that a nearly critically magne- 
tized dense core, ideally coupled to the magnetic field, is 
difficult to fragment into a few objects even if relatively 
large fluctuations are present initially. This conclusion 
is in good agreement with the conclusion reached by 
Hosking & Whitworth (2004), that a magnetically sub- 
critical cloud, which collapses after sufficient magnetic 
flux has been lost because of ambipolar diffusion, does 
not fragment. 

We conclude that, with large amplitude density fluc- 
tuations, a collapsing dense core fragments over a large 
range of initial conditions as long as it is not too strongly 
magnetized. The fragments are due to each initial seed 
collapsing individually. The question of the survival of the 
fragments should however be carefully investigated. We 
note that Price & Bate (2007) reach very similar conclu- 
sions. They also propose that in some cases, as when the 
magnetic field is perpendicular to the initial rotation axis, 
the magnetic field may help the fragment to survive. 

5. Discussion 

Since there is no question that a large fraction of stars are 
binaries, the efficient stabilization of the disk and the sup- 
pression of fragmentation even for modest field amplitudes 
may constitute a severe difficulty. We have no answer to 
this apparent conundrum but we discuss three possible 
mechanisms which may induce the fragmentation of the 
dense cores and the formation of binaries. We also sug- 
gest various observational tests which could discriminate 
between these scenarios. 

5.1. Large initial density fluctuations 

We start with the possibility of having sufficiently large 
fluctuations initially. As demonstrated in this paper and 
in Price & Bate (2007), fragmentation is possible in such a 
situation. The question is then how likely are large initial 
perturbations ? Observationally the question is not easy to 
address. It is nevertheless known that the dense cores usu- 
ally are not very uniform and that typically sonic velocity 
dispersion is observed. This is broadly compatible with 
the presence of density perturbations although maybe not 
as large as 50%. On the other hand, cores with a high 
aspect ratio are observed. If some of them are elongated, 
i.e. prolate objects, this may be equivalent to a substan- 
tial to = 2 perturbation. Indeed, simulations of quiescent 
cores having transonic turbulence, present significant fluc- 
tuations (Goodwin et al. 2004). It seems however difficult 
at this stage to infer quantitatively whether the perturba- 
tions generated by this weak turbulence are sufficient to 
produce multiple systems. 



10 



Hennebelle P., Teyssier R.: Magnetic processes in a collapsing dense core. II- Fragmentation 



Another possibility is that the collapse may be induced 
by an external agent such as supernova remnant or pro- 
tostellar jet. These are likely to generate large perturba- 
tions. Indeed, all simulations considering clouds evolving 
far from equilibrium (Bate et al. 2003, Ballesteros-Paredes 
et al. 2003, Hennebelle et al. 2006, Peretto et al. 2007) do 
find dense cores with initially strong perturbations. 

Future works should specify under which conditions 
such initial fluctuations are likely to be produced. This 
may depend on the physical properties of the molecular 
cloud in which the dense core is embedded. For exam- 
ple, the answer could be different in the Taurus molecular 
cloud in which star formation is relatively quiescent and 
in the Orion molecular cloud where star formation is more 
active. 

We stress that this mechanism constitutes a change 
of paradigm with respect to the standard hydrodynamical 
scenario (corresponding to Fig. 1). Indeed this mechanism 
implies that the fragmentation is the result of perturba- 
tions seeded at large scales rather than an intrinsic prop- 
erties of the collapsing dense core. 

An important prediction of this model is that the bi- 
nary should not be necessarily located in the equatorial 
plane of the core since the initial perturbations should 
likely be randomly distributed. 



value of Bg such that Bg/ \fAnp ~ C s , the ambipolar dif- 
fusion time becomes: 



Tad 



4nC 2 



(9) 



Following Shu et al. (1987), we write pi — C^fp 
where C = 3 x 10" 16 cm" 3 / 2 g 1 / 2 . With h ~ r/V2d, 
V e = y/M s G/r and B z = H z C 2 /VGr, we obtain 



Tad 



Af , H z 



(47r)3/ 2 C s V r V2d 



Thus, we obtain 



Tad 



= 3 X 



0.1M 



100AU H z 
r Vd 



(10) 



(11) 



For p = 20, H z ~ 0.4. From Fig. 2 of paper I, d ~ 10. 
Thus, the ambipolar diffusion time is comparable to the 
growth time of the toroidal component. It is therefore 
likely ambipolar diffusion will change the picture slightly. 
It may help clouds with low p to fragment more easily. 
We do not expect it however, to shift very significantly 
the value of p for which fragmentation is occurring. Since 
H z is proportional to l//i, T a d/T mag increases rapidly when 
P decreases. 



5.2. Ambipolar diffusion 

In the present calculations, perfect coupling between gas 
and magnetic field is assumed. The loss of flux due to am- 
bipolar diffusion could therefore possibly help to reduce 
the magnetic intensity. The ratio of the ambipolar diffu- 
sion time, r a d, and the freefall time, r//, has been com- 
puted by Shu et al. (1987). When the magnetic field is just 
sufficient to compensate gravity, they obtain T a d/Tff — 8, 
indicating that the ambipolar diffusion time is larger than 
the freefall time. When the magnetic field is weaker, the 
ambipolar time increases further. Therefore, it appears 
unlikely that magnetic field could be removed from the 
envelope of a supercritical collapsing core. 

However, the possibility remains that ambipolar diffu- 
sion could occur in the disk allowing the toroidal magnetic 
field to be diffused out. To investigate this issue, we es- 
timate the ratio of ambipolar diffusion time in the disk 
along the vertical direction and the growth time of the 
toroidal magnetic field given by Eq. ([2]). Here we write it 
as 



(7) 



B z V e 



The ambipolar diffusion time of the toroidal compo- 
nent along the vertical direction is: 



Tad 



IPPih 2 
Bl 



(8) 



where 7 = 3.5 x 10 13 cm 3 g 1 s 1 is the drag coefficient 
and pi the density of ions. Since we are interested in the 



5.3. Fragmentation during the second collapse 

Fragmentation during the second collapse has been inves- 
tigated by Bonnell & Bate (1994) with the prospect of 
exploring whether the formation of close binary systems 
could be possible. They conclude that the fragmentation 
is indeed possible but that the binaries have to accrete 
most of their mass since they are initially very small. 

On the other hand, Machida et al. (2007a) recently ex- 
plored the magnetized second collapse taking into account 
the large Ohmic dissipation which is predicted to occur in 
the first Larson core by various models (e.g. Nakano et al. 
2002). As a result, most of the magnetic flux is lost mak- 
ing the fragmentation easier to proceed. Indeed, the very 
recent study of Machida et al. (2007b) shows that frag- 
mentation during second collapse is not only possible but 
a promising mechanism. Interestingly, Banerjee & Pudritz 
(2006) also reports the formation of a very close binary in 
their second collapse calculations despite the ideal mhd 
assumption. It therefore appears possible that fragmenta- 
tion could occur during this phase possibly driven by rota- 
tion. In this case, the binary should gain sufficient angular 
momentum from the accretion to increase the separation 
between the two stars (see e.g. Goodwin et al. 2004). 

The prediction of this model is that the binaries should 
be in the equatorial plane since the orbital angular mo- 
mentum of the binary is due to the angular momentum 
of the accreting gas. Also the separation between the 2 
stars should increase with time implying that closer bi- 
naries should be observed on average into younger cores. 
Note that this mechanism could possibly work even if the 
magnetic field is initially very strong in the core. 
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In order to test each of these scenarios, high resolution 
observations would be necessary. In this respect, ALMA 
will certainly be a very powerful tool as demonstrated 
by the synthetic observations performed by Andre et al. 
(2007) using the present calculations. 

6. Conclusion 

In this paper we have studied the fragmentation of a col- 
lapsing magnetized molecular dense core by performing a 
set of numerical simulations. We interpret the results with 
the analysis performed in the companion paper in which 
we investigate the accretion and ejection processes that 
take place in the simulations. 

With our choice of initial conditions, and for pertur- 
bations of low amplitude, unmagnetized cores fragment, 
while magnetized cores having values of p, as large as 20 do 
not fragment. Based on an analytic estimate, we suggest 
that for cores having /x larger than ~ 15/a, the growth 
rate of the toroidal component is too slow to stabilize the 
disk. We stress that the suppression of fragmentation in 
this range of parameters, is not due to magnetic braking 
but to the rapid growth of the toroidal component of the 
magnetic field induced by the differential rotation within 
the disk. The Alfven speed associated to this toroidal com- 
ponent adds up to the sound speed of the disk and stabilize 
it. 

For values of /i smaller than 5, no big centrifugal 
disk forms because first the collapse occurs mainly along 
the field lines bringing less angular momentum and sec- 
ond magnetic braking removes angular momentum. This 
makes these cores even less prone to fragment. 

The situation is different if large amplitude perturba- 
tions are initially seeded. In this case, each perturbation 
develops independently even without rotation. Since the 
strong field amplification is primarily due to differential 
motions in the disk, the magnetic field is unable to sup- 
press fragmentation except if the core is almost critical, 
i.e. the field is initially very strong. The following evolu- 
tion of these fragments requires a careful treatment of the 
thermodynamics of the first Larson core. It is indeed likely 
that if the second collapse has not occurred by the time 
where these fragments approach each other, they are go- 
ing to merge. On the other hand, if the protostars already 
formed, the binary system is likely to survive. 

For dense cores having rotation and magnetic strength 
typical of values inferred from observations, we find that 
fragmentation is suppressed by the magnetic field if the 
initial density perturbations are too small. This consti- 
tutes a severe problem, since there is no question that 
a significant fraction of stars are binaries. In view of this, 
we discuss the likelihood of having sufficient perturbations 
within the cores initially, the impact of ambipolar diffu- 
sion and the possibility of fragmenting during the second 
collapse phase. We speculate that the first may depend on 
the physical characteristic of the molecular cloud in which 
the dense core is embedded whereas the last should be 



relatively independent of the large scales and could work 
even if the magnetic field is initially very strong. 
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